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We consider the generalization of the Kramers escape over a barrier 
problem to the case of a long chain molecule. The problem involves the 
motion of a chain molecule of segments across a region where the 
free energy per segment is higher, so that it has to cross a barrier. We 
consider the limit where the length of the molecule is much larger than 
the width of the barrier. The width is taken to be sufficiently wide that 
a coninuum description is applicable to even the portion over the barrier. 
We use the Rouse model and analyse the mechanism of crossing a barrier. 
There can be two dominant mechanisms. They are: end crossing and 
hairpin crossing. We find the free energy of activation for the hairpin 
crossing is two times that for end crossing. In both cases, the activation 
energy has a square root dependence on the temperature T, leading to 
a non-Arrhenius form for the rate. We also show that there is a special 
time dependent solution of the model, which corresponds to a kink in 
the chain, confined to the region of the barrier. The movement of the 
polymer from one side to the other is equivalent to the motion of the kink 
on the chain in the reverse direction. If there is no free energy difference 
between the two sides of the barrier, then the kink moves by diffusion 
and the time of crossing t cross ~ N'^/T^^'^. If there is a free energy 
difference, then the kink moves with a non-zero velocity from the lower 
free energy side to the other, leading to tcross ~ N/ ^/T. We also consider 
the translocation of hydrophilic polypeptides across hydrophobic pores, 
a process that is quite common in biological systems. Biological systems 
accompolish this by having a hydrophobic signal sequence at the end 
that goes in first. We find that for such a molecule, the transition state 
resembles a hook, and this is in agreement with presently accepted view 
in cell biology. 

I. INTRODUCTION 

The escape of a particle over a one dimensional barrier has been the subject of a large 
number of investigations. The problem, often referred to as the Kramers problem has 
been the subject of detailed reviews Kramers found solutions in the limit of weak 

friction and also in the limit of moderate to strong damping 0]. The intermediate regime 
has been an active area of investigation 0. The reason for this extensive activity is that this 
forms a model for a chemical reaction occuring in a condensed medium. Kramers problem 
for few degrees of freedom has also been the topic of study The quantum problem 

of escape/tunneling through a barrier too is of considerable interest. In the case where 
the system has an infinite number of degrees of freedom, this has been referred to as the 
decay of metastable vaccum, a problem that has attracted quite a bit of atttention in field 
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theory, cosmology and mesoscopic quantum phenomena In this paper, we consider a 

similar situation involving only classical physics. The trapped object has (— > oo) degrees 
of freedom, and is a polymer (a string). Though there are no quantum effects, the problem 
is similar, and some experiments are already available, that the results of the theory are 
easily verified. Further, the mathematics is considerably simpler than in the other cases, 
being equivalent to that of quantum mechanical tunneling of a single particle in a bistable 
potential. 

The way that the N degrees of freedom are connected (a chain or a string) leads to 
interesting new aspects to the problem that are not present in the case where there are 
only finite number of degrees of freedom. Also, the problem is of great interest in biology 
as many biological processes involve the translocation of a chain molecule from one side of 
a membrane to the other, through a pore in the membrane. The translocation of proteins 
from the cytosol into the endoplasmic reticulum, or into mitochondria or chloroplasts are 
processes of great interest and importance. Often, the proteins are hydrophilic and the pore 
in the membrane forms a hydrophobic region, through which it has to pass through [0-0, 
resulting in an increase in the free energy for the portion of the chain inside the pore. In 
infection by bacteriophages, conjugative DNA transfer etc, long chain DNA molecules snake 
through pores in membranes ||TT|JT^ . In all these cases, the chain molecule seems to get across 
the membrane rather easily, contrary to the expectation that one gets from the theoretical 



analysis available in the literature on the subject (see below). Bezerukov et. al. |jT5[ have 



studied the partitioning of polymer molecule into a nanoscale pore. Chipot and Pohoille [14 



have carried out a molecular dynamics simulation of a polypeptide, translocating through 
the interface between hexane and water. They found that the polypeptide (undecamer of 
poly-L-leucine) , initially placed in a random coil conformation on the aqueous side of the 
interface rapidly translocates to the interfacial region and then folds. In another interesting 



experiment, Han et. al. [jT3| observed the forced movement of long, double stranded DNA 
molecules through microfabricated channels which have regions that present an entropic 
barrier for the entry of the molecules. 

All these problems involve the passage of a long chain molecule, through a region in space, 
where the free energy per segment is higher, thus effectively presenting a barrier for the 
motion of the molecule. This problem forms the generalization that we refer to as the 
Kramers problem for a chain molecule. On the theoretical side, a variety of studies exist 
on this kind of problem. Muthukumar and Baumgartner [jl6| studied the movement of self 
avoiding polymer molecules between periodic cubic cavities seperated by bottlenecks. The 
passage through the bottleneck presents an entropic barrier to the motion, and they show 
that it leads to an exponetial slowing down of diffusion with the number of segments 



in the chain. Baumgartner and Skolnick [17] studied the movement of polymers through a 



membrane driven by an external bias and membrane asymmetry. Park and Sung P,p!8| , have 
studied the translocation through a pore. They analyze the passage through a pore on a fiat 
membrane, with only the effects of entropy included. The resultant entropic barrier is rather 
broad, its width being proportional to A^. Consequently, they consider the translocation 
process as being equivalent to the motion of the center of mass of the molecule. Using the 
result of the Rouse model that the diffusion coefficient of the center of mass is proportional 
to 1/A^, they effectively reduce the problem to the barrier crossing of single particle having a 
diffusion coefficient proportional to 1/A^. As the translocation involve motion of segments 
across the pore, the time taken to cross, tcross scales as A^^. They also show that in cases 
where there is adsorption on the trans side, translocation is favored and then tcross scales 
as A^^. In a very recent paper. Park and Sung [PTI] have given a detailed investigation of 
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the dynamics of a polymer surmounting a potential barrier. They use multidimensional 
barrier crossing theory to study the motion of a chain molecule over a barrier, in the limit 
where the width of the barrier is much larger than the lateral dimension of the molecule. 
In an interesting recent paper, Lubensky and Nelson study a case where they assume 
the interaction of the segments of the polymer with the pore to be strong. They argue 
that effectively, the dynamics of the portion of the chain inside the pore is the one that is 
important and this, they show, can give rise to tcross proportional to A^. Again, they assume 
diffusive dynamics. In a recent paper, we have suggested ||23| a kink mechanism for the 
motion of the chain across a barrier and it is our aim to give details of this mechanism in 
this paper. 

We consider a polymer undergoing activated crossing over a barrier. This can form a 
model for a polymer going through a pore too, as the pore can cause an increase in the free 
energy of the segments inside it, as they would interact with the walls of the pore. The 
width w of the barrier is assumed to be much larger than the Kuhn length / of the polymer, 
but small in comparison with the total length Nl of the polymer. That is, I << w << Nl 
. For example, in eukaryotic cells 0, the length of the nuclear pore is about 100 A, while 
the Kuhn length for a single stranded DNA is perhaps around 15 A |2^. Therefore, one is 



justified in using a continuum approach to the dynamics of the long chain. (It is possible 
to retain the discrete approach, and develop the ideas based on them, but this is more 
involved mathematically). Our approach is the following: We describe the motion of the 
polymer using the Rouse model. The force that the barrier exerts on the chain appears as 
an additional, non-linear term in the model. We refer to this as the non-linear Rouse model. 
The non-linear term causes a distortion of the portion of the chain inside the barrier, which 
we refer to as the kink. Movement of the chain across the barrier is equivalent to the motion 
of the kink in the reverse direction. The kink is actually a special solution of the non-linear 
Rouse model, arising because of the non-linearity. In the presence of a free energy difference 
between the two sides, the kink moves with a definite velocity and hence the polymer would 
cross the barrier with tcross proportional to A^. Traditionally, the non-linear models that 
one studies (for example, the 0^or the sine-Gordon model P^^) have potentials that are 
translationally invariant, and hence the kink can migrate freely in space. In comparison, in 
our problem, the non-linear term is fixed in position space and hence the kink too is fixed in 
space. However, the chain molecule (modelled as a string [|l9l) can move in space and hence 
the kink migrates, not in space, but on the chain. As far as we know, such a suggestion has 
never been made in the past and we believe that this is a very useful idea in understanding 
polymer translocation. 

In general, the polymer can escape by essentially two mechanisms. The first, which we 
refer to as end crossing, involves the passage of one end of the polymer over the barrier, 
by thermal activation. This leads to the formation of the kink, which is then driven the 
free energy difference between the two sides of the barrier. The second is by the escape 
of any portion of the polymer over the barrier, in the form of a hairpin. The hairpin is 
a kink-antikink pair. For a flexible polymer, the hairpin crossing has twice the activation 
energy for end escape and hence one expects it to be less probable. However, as it can take 
place anywhere on the chain, the frequency factor for it is proportional to A^ and hence for 
a sufficiently long chain, this can become the dominant mechanism for the escape. Hairpin 
crossing leads to the formation of a kink-antikink pair. The pair moves apart on the chain, 
driven by the free energy gain and hence the time of crossing is still proportional to A^, 
though one expects that it is roughly half the time of crossing in the end crossing case. 
In principle, in addition to these, it is possible for more than one hairpin to be formed. 



3 



However, it is obvious that in passage through a pore, unless the pore is rather wide, only 
the end-crossing mechanism would operate. 

In all our mathematical development, we use the one dimensional version of the Rouse 
model. This is no limitation, if one is concerned with translocation across the interface 
between two immiscible liquids or the experiments of Han et. al. |[T^, which involve motion 
in a channel, whose width is large in comparison with the size of the molecule. On the 
other hand, if one is interested in translocation through a pore, strictly speaking, one has 
to consider the full three dimensional nature of the problem, which at present seems rather 
involved. However, we believe that the one dimensional model captures the essential physics 
of the problem. Our analysis should also be useful in situtations where the whole of the 
polymer is in a pore, so that the dynamics may be taken to be one dimensional, with the 
chain trying to cross a region of high free energy. 



II. THE MODEL 
A. The Free Energy Landscape 

The considerations in this section are quite general and do not depend on the model that 
one uses to describe the polymer dynamics but involves the assumption that the polymer is 
flexible over a length scale comparable to the width of the barrier. We start by considering 
the free energy landscape for the crossing of the barrier. The barrier and the polymer 
stretched across it are shown in the figure The polymer has initially all its units on the cis 
side, where its free energy per segment is taken to be zero. So the initial state, corresponds 
to a free energy of zero in the free energy hypersurface shown in figure 0. In crossing over 
to the trans side, it has to go over a barrier, as in the figure |1|. The transition state for the 
crossing can be easily found, from physical considerations. It is the state shown in figure 
^. In it, the configuration of the polymer is such that the end of the polymer on the trans 
side is located exactly at the point on the trans side at which its free energy per segment 
is zero, with the other end on the cis side, and the chain is such that the free energy of the 
whole chain is a minimum. This is so because if one moves the end either in the forward or 
in the backward direction (and the rest of the chain adjusted so that the free energy of the 
chain as a whole minimum), then the total free energy of the system would decrease. Hence 
in the free energy hypersurface figure |^, the configuration shown in the figure ^ corresponds 
to the maximum (i.e. transition state). Once the system has crossed the transition state, 
the chain is stretched across the barrier. The path of steepest descent then corresponds to 
moving segments from the cis side to the trans side, with out changing the configuration of 
the polymer in the barrier region. As there is a free energy difference AV between the two 
sides, this would lead to a lowering of the free energy by AV per segment, and this leads to a 
path on the free energy surface with a constant slope, and of width W proportional to (see 
figure 3). Such a landscape implies that the translocation process would involve two steps. 
First step is going through the transition state by the overcoming of the activation barrier. 
Once the system has done this, it encounters a rather wide region of width proportional 
to the length of the chain. Crossing this region is the second step. As this region has a 
constant slope, the motion in this region is driven and it is similar to that of a Brownian 
particle subject to a constant force. Such a particle would take a time tcross, proportional to 
A^ to cross this region. 

Till now, we considered the case of end-crossing. The scenario for hairpin crossing is 
similar. However, the activation energy is higher for hairpin crossing. In hairpin crossing. 
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the transition state is equivalent to the one end crossing, repeated two times. Hence the 
activation energy for the process is two times larger. Once a hairpin crossing occur, a kink- 
antikink pair is formed and the kink and the anti-kink separate rapidly, due to the driving 
force of the free energy gain. Then futher crossing occurs by the movement of these two on 
the chain, which again leads to a time of crossing proportional to A^. 

In the following we make all these considerations quantitative, using the Rouse model to 
describe the dynamics of the chain. 



B. The Dynamics 



We consider the continuum limit of the Rouse model, discussed in detail by Doi and 
Edwards The chain is approximated as a string, with segments (beads) labelled by 

their position n along the chain, n is taken to be a continuous variable, having values 
ranging from to A^. The position of the ra*^ segment in space is denoted by R{n,t), where 
t is time. In the Rouse model, this position undergoes overdamped Brownian motion and 
its time development is described by the equation 



c 



dR{n,t) d'^R{n,t) 



dt 



m- 



V'{R{n,t)) + f{n,t). 



(1) 



In the above, C is a friction coefficient for the n*^ segment. The term al"'*'' comes from 
the fact that stretching the chain can lower its entropy and hence increase its free energy. 
Consequently, the parameter m = SksT/P (see Doi and Edwards equation (4.5). They 
use the symbol k for the quantity that we call m) . As the ends of the string are free, the 
boundary conditions to be satisfied are / 9R{n,t) j _ f 8R(n,t) j _ V(R) is the free 

J I dn )n=o I }n=N ^ ' 

energy of a segment of chain, located at the position R. We assume that V[R) leads to a 
barrier located near i? = 0. f{n,t) are random forces acting on the n'^^ segment and have the 
correlation function (/(n, t)/(ni, ti)) = 2(kBTS{n — ni)6{t — ti)(see [|19[, equation (4.12)). 
The deterministic part of the equation (|l]), which will play a key role in our analysis, is 
obtained by neglecting the random noise term in (|l|). It is: 



c 



dR{n,t) d'^R{n,t) 



dt 



m- 



V'{R{n,t)) 



This may also be written as: 



c 



dR{n,t) 5E[R{n,t)] 



dt 6R{n, t) 

where E[R{n,t)] is the free energy functional for the chain given by: 



(2) 



(3) 



E[R{n,t)] 



N 



dn 



m ( dR{n, t) 
~2 \ 



+ V{R{n,t)) 



(4) 



C. The form of the barrier 



The chain is assumed to be subject to a biased double well potential (BDW), of the form 
shown in the figure 1. The two minima are at — ao and ai, with < ai. There is assumed to 
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be a maximum at i? = 0. Further, we take V{—ao) = 0. All these conditions can be satisfied 
if one takes V'{R) = 2 k R{R + oq) {R — cti). Here, /c is a constant and will determine the 
height of the barrier. Integrating this and using V{—ao) = 0, we get 

k 

V{R) = -{R + aof{3R^ - 2Rao - ARai + + 2aoai) (5) 
6 

The barrier height for the forward crossing is Vf = V{0) — V(— Oq) = ^kaQ^ (oq + 2ai) and 
for the reverse process, it is VJ, = V{0) — V{ai) = ^kai^ (2ao + cti). On crossing the barrier, a 
unit of the polymer lowers its free energy by AV = V"(ai)— V(— oq) = ^k (ag — oi) (o-o + ^i)^. 
The form of the potential is shown in the figure |^. 



D. The Activation Free Energy for End and Hairpin Crossings 



In this section, we consider the first step and calculate the activation free energy for both 
end and hairpin crossing. Activation free energy can be obtained from the free energy func- 
tional of equation . This free energy functional implies that at equilibrium, the probability 



distribution functional is exp 

polymer 
equation 



Idnnm{^) +V{R{n)) 



The configurations of the 



polymer which makes free energy a minimum are found from "^^^j^^^ = 0, which leads to the 

- V'W (6) 

Notice that this is just a Newton's equation for a particle (ficticious, ofcourse) of mass m 
moving in a potential —V{R). This equation has four solutions that are of interest to us. 
The first two are: (1) R{n) = —aQ, (2) R{n) = ai which are the minima of the free energy. 
The first solution is the initial state, where the polymer is trapped in the vicinity of — oq. 
The second is the most stable minimum, at R{n) = ai. In addition to these, there are two 
more solutions which are of interest to us. These are n dependent and correspond to end 
and hairpin crossings. 



1. End Crossing 

As we are interested in the case where the polymer is very long, we can imagine n to 
vary from — oo to and find a saddle point in the free energy surface by searching for 
a solution satisfying R{—oo) = —ao and the other end of the polymer to be at a point 
with R > -Rmax, where -Rmax is the point where V{R) has its maximum value. For the 

Newton's equation (||) the conserved energy is Ec = \m —V{R{n)). For the extremum 
path, Ec = 0. Thus, the particle starts at R{—oo) = — ao with the velocity zero (this 
follows from the boundary conditons of the Rouse model) and ends up at i?/ at the "time" 
n = 0. Here -R/(> -Rmax), is the point such that V{Rf) = 0, again with the velocity zero. 
Further, free energy of this configuration is activation free energy for end crossing. As for 

this configuration, = V{R{n)), we find the activation free energy to be given by 

Ea,end= f^' ^2mV{R)dR. (7) 

J-ao 

The end crossing is illustrated in figure |^. 
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2. Hairpin Crossing 



If one imagines n to vary in the range (—00, 00) a second saddle point may be found by 
taking R{—oo) = — ao and R{oo) = — ao, so that the Newtonian particle starts at — ao, 
makes a round trip in the inverted potential —V{R) and gets back to its starting point. This 
obviously has an activation energy 



Ea,hp = 2 pmV{R)dR = 2Ea,end 

J— an 



(8) 



Thus the activation energy is exactly two times for end crossing pO[- The hairpin crossing 
is shown in figure |^ 



3. The Temperature dependence 

As the parameter m is proportional to the temperature ( = SksT/P ), we arrive at the 
general conclusion that both the activation energies Ea^end and Ea^hp are proportional to \/T. 
For our model potential of equation we find Rf = ao(7 — 1/7^ — 7) where 7 = (1 + 2^)| 
and 



E. 



a, end 



'mkal 



(3 7' + 1)^1 + 37 - 3 7 (7' - 1) In (v/7(7 - 1)/ (l + 7 - / 



1 + 37 



(9) 



The Boltzmann factor e ''b'^ for the crossing of one end of the polymer over the barrier thus 
has the form e-constant/^/^_ Pui-^j^er, we find that it is independent of for large N. 



III. THE RATE OF CROSSING 



A. Hairpin Crossing 



We now calculate the rate of crossing in the two cases. We first consider the hairpin 
crossing, as this has connections with material available in the literature The methods 
that we use are quite well known in the soliton literature P3 and hence we give just enough 
details to make the approach clear. The Rouse model in the equation ([l|) leads to the 
functional Fokker Planck equation 



dP 
'di 



= — dn 



N 



6 



C Jo 6R{n) 



knT 



6P 



6R{n) 6R{n 



(10) 



for the probability distribution functional P. This equation implies that the flux associated 
with the co-ordinate R{n) is |^| 



j{R{n)) = - 



c 



6P 



SR{n) 5R{n 



'IV. 



We now consider the initial, metastable state. As the rate of escape is small, we can assume 
the probability distribution to be the equilibrium one, which is 
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P=^e^^{-E[R{n)]/kBT} (12) 

To determine Zq we use the condition / D[R{n)] P = 1, where / D[R{n)] stands for functional 
integration. It is convenient to introduce the normal co-ordinates for small amplitude motion 
around the metastable minimum and do the functional integration using them. For this, we 
expand E[R{n)] around the metastable minimum, by putting R{n) = — Oq + 6R{n), and 
expanding as a functional Taylor series in 6R{n) and keeping terms up to second order in 
6R{n). Then 

E[R{n)] = -m 6R{n) (^-— + u^j SR{n) (13) 

We have defined Uq by putting muQ = ^ YSi'^ . The normal (Rouse) modes are just 

the eigenf unctions ipki^) of the operator H"^^ = (— ^ + i-Uq^, having the eigenvalue Ek and 
satisfying the Rouse boundary conditions = at the two ends of the string. (The 

superscript "ms" in i^'^^stands for metastable). Now we can expand 5R{n) as 5R{n) = 
J^kCk'ipkin) so that the expression for energy (^) becomes 

E[R{n)] = ^mJ2ekcl (14) 

k 

We now do the functional integration using the variables Ck- Then the normalization 
condition J D[R(n)] P = 1 becomes J- 11 J dckexp —^mPekcl = 1. This leads to 

k 

/ 9 \l/2 

= II I a ■ Now we consider the vicinity of a saddle point, where the probabil- 

ity distribution deviates from the equilibrium one. We first consider the saddle point which 
corresponds to hairpin crossing. The potential of the equation (^ is rather difficult to han- 
dle as we have not been able to obtain analytic solutions to the Newton's equation (^. 
In determining the crossing of the barrier, the key role is played by the quantities ujq and 
the height of the barrier for crossing in the forward direction Vj. The quantities that we 
calculate in this section have no dependence of the behavior of the potential near the stable 
minimum. So, instead of using the quartic potential of the equation (^, we use the simpler 
cubic potential of equation (^). This has no stable minimum (corresponding to the final 
state), but that does not matter, because the quantities that we calculate do not depend on 
its existence. Thus we use the potential: 

K(ii) = V'„(^)'(l-^) (15) 

where we adjust Vq and Rq to reproduce the values for uq and the barrier height V}. Solving 
the equation @ for this potential, in the limit of an infinitely long chain extending from 
n = —00 to +00, the saddle point that corresponds to hairpin crossing is easily found to be 
given by the equation 

Rhpin) = -ao + Ro jsec h ^y'^'^j | (16) 

In fact one has a continuous family of solutions of the form Rhp{n — no), where no € {—00, 00) 
is arbitrary and determines center of the kink-antikink pair. Now expanding the energy 
E[R{n)] about this saddle, by writing R{n) = Rhpiji — no)+ 5R{n) we get 
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E[R{n)] = Ea,hp + J dn6R{n) + cug |l - 3 sec /i^ (luoI^ - no)/2)| 



SR{n) (17) 



For the potential of equation ( pTS]) -Ea,/ip = (8-Ro/15) ^/2rnVo. The normal modes for fluc- 
tuations around the saddle are determined by the eigenfunctions of the operator = 
— ^ + cJq |l — Ssech^ {iOo{n — rio)/2)| . (t is used to denote the saddle point). The eigen- 
functions are: 
and el 



(aj the discrete states tpl, ipl and ipl having the eigenvalues el = —bu^/A, 
3^0/4 and (b) the continuum of eigenstates with eigenvalues of the form 
el = uJq + fc2(more details are given in appendix A). We denote the eigenfunctions by ipl- 
The existence of the eigenvalue £1 = comes from the freedom of the kink-antikink pair to 
have its center anywhere on the chain (the hairpin can be formed anywhere). In the follow- 
ing, would stand for summation over all the eigenstates, including both the discrete and 

k 

continuum states while a symbol like J2 means that the bound state ipl is to be excluded 

k^l 



4 



el 



from the sum. Now writing 6R{n) 



E cI^IjI, we get 



E[R{n)] = E,,,, + ^mJ2ei{cl)' 
We write the probability density near the saddle as 

P 



OiclcU) f E[R{n)] 
exp <. 



where 9{cl,Ci...), is a function that must approach unity in the vicinity of the metastable 



minimum. Near the saddle, one can calculate the flux jl in the direction of c] 



i 



\^r^dP_ ^ dE[R{n)] ^ 
dc\ dc\ 



Using the equations (|T^ ) and ([T^ ) we get 

kBTde{cic\...) 



del 



exp 



(19) 



In a steady state, there is flux only in the unstable direction. That is, only Jq is non-zero. 
This means that 6 can depend only on Cq, which impies that jg must have the form 



il = A exp 



ksT \ 2 



-m 



E 

fe>i 



ei c 



(20) 



where A is a constant, to be determined. Using the equation (pO|) in ([19D we get 



del 



—A exp 



2kBT 

one to get A = 



Cq) I. The fact that 0{cl) must approach unity as Cq — > —00, enables 



ti \ 1/2 



2nkBT 



. Hence 6{q 



ti \ 1/2 



2-wkBT 



/^Trfzexp{-2^m 4 z'^]- Now 
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the net flux crossing the barrier is found by integrating jg over all directions other than 
Cq . The integrals over all c|., except cf is straightforward. As = 0, / del needs special 
handling. The integral, as is well-known, is performed by converting it to an integral over the 

kink-antikink position, uq. That is, / dc\ = aj dno, where = J^^^dn (^-^||^) = 
Hence the rate becomes 



,,hp 



m 



hp 



1/2 



n 



2-KkBT 



1/2 



-t 



Ea. 



hp 



m 



1/2 



NeK^{-Ea,hp/kBT) 



(21) 



The notation H is used to indicate product over all eigenvalues of H^, except the first two. 
On using the expression for Zq, 



khp — 



where Ihp 




kBT / m \3/2 

1/2 



4 


Ea,hp 




=-2 


m 



1/2 



Nexp{-Ea,hp/kBT) 



(22) 



This infinite product is evaluated in the appendix B and is found 



> fe>2 

to be Ihp — Y^o- This leads to 



^hp 



47rC \27rkBT 



exp {-Ea,hp/kBT) 



(23) 



B. End Crossing 

In this case, the analysis is similar to the above. The operator is the same as earlier. 
However, there is an interesting difference. In the hairpin case, the boundary conditions on 

( ^ = 0, at the two ends) were at n = ±oo, while in this case, they are at n = and 
at n = oo (i.e. the boundary value problem is now on the half- line). Due to this, one has to 
rule out the odd that exists in the hairpin case as they do not satisfy the Rouse boundary 

condition-^ = at n = 0. So we consider only the even solutions. Thus the eigenvalue at 
zero is ruled out (which is quite alright as end crossing can occur only at the end and not 
anywhere else, but we will put in additional factor of 2 as it can occur at the two ends). 
The discrete spectrum now has only the eigenvalues s\ = —buj^/A, and £\ = 3a;Q/4. The 
expression for the rate is 



_ kBT I rn 



C I 2T:kBT 



1/2 

lend exp (-£^a,end/ ^S^) (24) 



. 1/2 



where lend = ■^"^ A ,2 . In this product, there are — 1 terms in the numerator and N 



terms in the denominator. One of the N — 1 terms is the bound state with an eigenvalue 
£2 = 3oUq/A. Separating this out from the product, one can write lend = (s^^l^U^) ^^nd, 
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rr \ 1/2 

k 
k>2 

in the Appendix B. The result is 



where lend = ( 1 ■ evaluation of this product involves some subtelity and is done 



5mC(jQ / r-i II m\ 

kend = — 7=-— exp {-Ea,end/kBT) 



Accounting for the existence of two ends leads to 



ktwo-ends — ^^P Ea,end/ ^bT) (25) 



IV. THE KINK AND ITS MOTION 
A. The kink solution and its velocity 

Having overcome the activation barrier, how much time would the polymer take to cross it? 
We denote this time by t cross- To calculate this, we first look at the mathematical solutions 
of the deterministic equation (|). The simplest solutions of this equation are: R{n,t) = — oq 
and R{n,t) = ai. These correspond to the polymer being on either side of the barrier and 
these are just mean values of the position on the two sides. Thermal noise makes R{n, t) 
fluctuate about the mean position which may be analyzed using the normal co-ordinates for 
fluctuations about this mean position. Each normal mode obeys a Langevin equation similar 
to that for a harmonic oscillator, executing Brownian motion. In addition to these two time 
independent solutions, the above equation has a time dependent solution (a kink) too, which 
corresponds to the polymer crossing the barrier. We analyze the dynamics of the chain, with 
the kink in it, using the normal modes for fluctuations about this kink configuration. Our 
analysis makes use of the techniques that have been used to study the diffusion of solitons 



As is usual in the theory of non-linear wave equations, a kink solution moving with a 
velocity v may be found using the ansatz R{n,t) = Rs{t) where t = n — vt [El]. Then the 
equation reduces to 

If one imagines r as time, then this too is a simple Newtonian equation for the motion of 
particle of mass m, moving in the upside down potential — V{R). However, in this case, there 
is a frictional term too, and vC,/m is the coefficient of friction. This term makes it possible 
for us to find a solution for quite general forms of potential, with V'{R) — as i? — > ±cxd. 
For the potential of the equation (^, we can easily find a solution of this equation, obeying 
the conditions Rs{t) = — oq for r — >• — oo and Rs{t) = ai for r — > oo. The solution is 

TUJ (ao+ai) \ ( -I I ™(ao+ai)\ 



R^{r) = -ao + e^" ^'^"+'^^^ ai 1 + e™ ^""+'^1^ , (27) 



with u = ^k/m . The solution exists only if the velocity v = (ag — Oi). This solution 
is a kink, occurring in the portion of the chain inside the barrier. We shall refer to the point 



11 



with r = as the center of the kink. (Actually one has a one-parameter family of solutions 
of the form Rs{t + Tq), where tq is any arbitrary contant). As t = n — vt, the center of 
the kink moves with a constant velocity v. Note that this velocity depends on the shape 
of the barrier. Thus for our model potential, if oq < ai, then Vf < Vb, and this velocity is 
negative. This implies that the kink is moving in the negative direction, which corresponds 
to the chain moving in the positive direction. That is, the chain moves to the lower free 
energy region, with this velocity. If the barrier is symmetric, then CLq — CLi 

( Vf = Vb) the 

velocity of the kink is zero. 



B. Fluctuations about the kink 

We now analyze the effect of the noise term present in the equation ([l|). The center of the 
kink can be anywhere on the chain - which means that the kink is free to move on the chain. 
Actually, as the position of the kink is fixed in space, this means that the polymer is moving 
across the barrier. The kink would also execute Brownian motion, due to the noise term. 
The motion of the kink caused by the noise terms is a well studied problem in the literature 



2^ and one can make use of these methods. Following Tnstanton methods' of field theory 



1 2^, we write 



R{n, t) = Rs{n - a{t)) + ^ Xp{t)(f)p{n - a{t),t) (28) 

p=i 

We have allowed for the motion of the kink by taking the kink center to be at a{t), where 
a{t) is a random function of time which is to be determined, (pp are a set of functions (the 
Rouse modes) below and Xp{t) are expansion coefficients. This may be put into the equation 
(|^) to derive an equation of motion for a{t). Neglecting kink-phonon scattering leads to 



a{t)=v + Ut)/C (29) 
where we define iJoi^n) by dnRsin) = Cipliji) with 

vC,n/m 



C = (dMn 

and 



e 



dnRs[n)) = —TTUJ csc(2 7r ) (ai — ao) aoOi- (30) 

'3 ttQ + ai 



1 /■N/2 ^ 

Ut) = 7 / dnroin)e<''/^''^^f{n + a{t),t). (31) 

C J-N/2 

^o{t) is a random function of time, having the correlation function 

mMh)) = Kt - h){2kBT/0 f , rfne™^/"^ [Mnf (32) 



lover the kink 

.For the potential given by the equation (^ one gets 

(^o(^)^o(^i)) = S{t- ti)kBT/ (2Cao ai) sec(2 tt "° ) (3 ai - ao) (3 ao - ai) . (33) 

The equations (|29|) and (|30D imply that the kink position a{t) executes Brownian motion 
with drift. As v is negative, the drift is in the negative direction. 
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C. The crossing time tcross 

For the polymer to cross the barrier, the kink has to go in the reverse direction, by a 
distance equal to A^. As the equation ( ]29|) is just that for a particle executing Brownian 
motion with drift, we can estimate the time of crossing as a first passage time. As the kink 
starts at one end, we take the initial position of the particle, a to be A^ and calculate the 
average time required for it to attain the value 0, which would correspond to the polymer 
crossing the barrier fully. Writing the diffusion equation for the survival probability P{a,t) 
for a particle starting at a = A^ at the time t = and being absorbed at a = 0, we get 



dP{a, t) 
di 



d^P{a,t) dP{a,t) 

L) V- 



da? 



da 



(34) 



Here, the diffusion coefficient 



D 



dti / dt2 {Uh)Uh)) 



SksT rm(3ai — a(j) (3ao — cti) ai — ao, 

a r \ T ^tan(27r^ — , 



(35) 



The equation (|^) is to be solved, subject to the initial condition P{a, 0) = 6{a — N) and 
with absorbing boundary condition at a = (i.e. P(0,t) = 0) and P(oo,t) = 0. It is easy 
to solve the above equation in the Laplace domain. The result for the Laplace transform 
P{a,s) = dtP{a,t) exp{—st) is: 



P{a,s) 



1 



(g-iV) v-\'4D 3 + \a-N\ (a-N)v-\/ AD s + v'^ \a\~\/ A D a + v'^ N 

e 2 U — g 2D 



(36) 



The Laplace transform of the survival probability is given by P{s) 
found to be 



P{s) 



-N v-yj AD s+v^ N 

1 — e 2D 



X!^ daP{a, s) and is 



(37) 



The average crossing time is given by tcross = Limit s->oP{s) = N/{—v) , if v < 0. As f is 
proportional a/ mfc, assuming V{R) to be temperature independent we find tcross ~ N/\/T. 
This is a general conclusion, independent of the model that we assume for the potential. If 
the barrier is symmetric, the kink moves with an average velocity v = Taking the f — >■ 
limit of -P(s), we get 



Pis) 



ill 

s 



so that the survival probability becomes 



P{t) = Erf( 



e 



N 



2VDt 



{31 



(39) 



This expression for the survival probability implies that the average time that the particle 
survives is tcross ~ N'^/D. For the symmetric barrier, the value of D may be obtained by 
taking the limit Oi ^ Oq, and one finds D = jf^^J^ and thus tcross ~ A^^/T^/^. 
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In their analysis, Park and Sung considered the passage of a polymer through a pore 
for which the barrier is entropic in origin. Consequently it is very broad, the width being 
of the order of N. Hence they consider the movement as effectively that of the center of 
mass of the polymer which diffuses with a coefficient proportional to As the center 

of mass has to cover a distance N, the time that it takes is proportional to A^^. If there 
is a free energy difference driving the chain from one side to the other, then the time is 
proportional to N"^. In comparison, we take the barrier to be extrinsic in origin and assume 
its width to be small in comparison with the length of the chain. The crossing occurs by the 
motion of the kink, which is a localized non-linear object in the chain whose width is of the 
same order as that of the barrier. As the polymer is intially subject to a potential well, the 
entropic contribution to the barrier that Park and Sung |jl8| consider does not exist in our 
case. Such a potential is realistic, in cases where the polymer is subjected to a driving force 
(for example an electric field). As the kink is a localized object, its diffusion coefficient has 
no A^ dependence and our results are different from those of Park and Sung . In the case 
where there is no free energy difference, our crossing time is proportional to A^^(in contrast 
to A^^ of Park and Sung) , while if there is a free energy difference, our crossing time is 



proportional to A^ (in contrast to A^ of Park and Sung). In a very recent paper |21|, Park 



and Sung have considered the Rouse dynamics of a short polymer surmounting a barrier. 
The size of the polymer is assumed to be small in comparison with the width of the potential 
barrier. Consequently, the transition state has almost all the units at the top of the barrier, 
leading to the prediction that the activation energy is proportional to A^. This leads to a 
crossing probability that decreases exponentially with A^. In comparison, as found in section 
11 D|, the free energy of activation does not depend on the length of the chain. Hence, the 



mechanism is the favoured one for long chains. 



D. The net rate 



As the actual crossing is a two step process, with activation as the first step and kink 
motion as the second step, the net rate of the two has to be a harmonic mean of the two 
rates. For a very long chain, the motion of the kink has to become rate determining. In 
the case of translocation of biological macromolecules, considered in section there does 
not seem to be any free energy of activation and then the rate is determined by tcross alone. 
Recently, the motion of long chains in microfabricated channels have been investigated by 
Han et al In contrast to the situtation for a pore, there is an additional direction is 



available for the molecule to form a hairpin, viz. perpendicular to the direction of movement 
of the molecule. Consequently, in overcoming the barrier, both end crossing and hairpin 
crossing can occur (see figures ^ and Experimental results show that the longer molecule 
crosses the barrier faster. This means that the A^-dependence of khp causes the hairpin 
crossing to be the dominant mechanism of crossing in these experiments. 



V. HOW DO BIOLOGICAL SYSTEMS LOWER THE ACTIVATION ENERGY? 

If there was a high activation energy (>> ksT) for the translocation, the process would 
be unlikely and hence, biological systems would not be able to function, if they depended 
crucially on such transfers. As translocation seem to be very efficient in biological systems, 
one needs to look at the mechanism that evolution has designed to reduce the barrier. The 
destination (referred to as sorting) of a biological long chain molecule is determined by 
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a sequence of units at the begining of the chain, referred to as the signal sequence. For 
example, proteins destined to the endoplasmic reticulum possess an amino-terminal signal 
sequence, while those destined to remain in the cytosol do not have this. If one attaches 
this sequence to a cytosolic protein, then the protein is found to end up in the endoplasmic 
reticulum (see reference figure 14.6). The way the sequence works is simple. If the pore 
is hydrophobic and the chain hydrophilic, then the signal sequence is hydrophobic, so that 
the signal sequence has a low free energy inside the pore. 

We qualitatively analyze this type of problem in the following, using the Rouse model. The 
way to model the situation would be to have a potential that is dependent upon the segment 
number n in the chain. Hence, in the equations of the Rouse model the potential term would 
have an explicit dependence on n. Let us denote the length of the signal sequence by s. The 
simplest model would be to have a potential which is attractive, for < n < s and which 
has the shape of a barrier for s < n < N. The transition state is determined by the the 
Newton-like equation 

(PR 

m-— = V^^Jn,R), (40) 

with n playing the role of time (in the following we shall refer to n as the time for the motion 
of this ficticious particle). We take the potential to be such that 



Vnew{n, R) = -V{R) if < n < s and 



= V{R) if s<n<N 

This corresponds to a particle moving in a time dependent potential, which switches from 
being repulsive to attractive at the time s. The shape of this time dependent potential is 
shown in the figure ||. The boundary conditions |^^^} ^ = {^^^} ~ ^ imply that 
the particle has to start and end with zero velocity. Let us imagine that the particle starts 
at the point Rq (see figure ||). As the potential that it feels up to the time s is repulsive, it 
follows the path indicated by the dashed line in the figure, and the conservation of energy 

may be written as: |m + ViR) = ViRo)- Let it reach the point Rs after a time 

s. At this time, the potential is switched from V{R) to —V{R). From this time on, the 
equation of conservation of energy would be: 




dR{n)y 
dn I 



V{R) = V{R^)-2V{Rs). (41) 



This is the equation of motion of the particle for s < n < N. We are interested in oo 
limit and we have to satisfy the boundary condition |^^^| ^ = at the end of the chain. 
In the particle picture, this is equivalent to the condition that the total energy of the particle 
obeying the equation ^ must be zero. This implies that V{Ro) = 2V{Rs). For a given s, 
this uniquely fixes the values of the two variables Rq and Rg. 

The net transition state is shaped like a hook and the hydrophobic part of the chain is 
completely in the short arm of the hook (see figure |A 11 2| ) . A configuration like the one in the 
figure [ly where the whole of the hook is formed by the hydrophobic part is not a transition 
state. The transition state in figure |A 11 2|, though it seems likely to occur in crossing between 
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liquid liquid interfaces, it seems rather difficult to form in the case of passage through a pore 
as there are two difficulties: (1) the chain has to bend to form the hook (2) the pore has 
to be wide enough to accommodate the two strands of the hook simultaneously. Inspite of 
these, nature does seem to use this as an inspection of the figure 14-14 of reference shows. 



VI. CONCLUSIONS 



We have considered the generalization of the Kramers escape over a barrier problem to 
the case of a long chain molecule. It involves the motion of chain molecule of segments 
across a region where the free energy per segment is higher, so that it has to cross a barrier. 
We consider the limit where the width of the barrier w is large in comparison with the 
Kuhn length Z, but small in comparison with the total length Nl of the molecule. The limit 
where Nl « w has been considered in a recent paper by Park and Sung [^]. We use 
the Rouse model and find there are two possible mechanism that can be important - end 
crossing and hairpin crossing. We calculate the free energy of activation for both and show 
that both have a square root dependence on the temperature T, leading to a non-Arrhenius 
form for the rate. We also find that the activation energy for hairpin crossing is two times 
the activation energy for end-crossing. Inspite of this, for long enough chains, where the 
geometry of the systems permits, hairpin formation can be the dominant mode of escape as 
seen in the experiments of Han et. al. [1^. The width of the barrier in these experiments 



is rather large in comparison with the length of the polyme so that the kink mechanism of 
crossing seems to be unlikely in this case. 

While in the short chain limit Park and Sung find the activation energy to be linearly 
dependent on A^, we find that for long chains, the activation energy is independent of A^. We 
also show that there is a special time dependent solution of the model, which corresponds to 
a kink in the chain, confined to the region of the barrier. In usual non-linear problems with 
a kink solution, the problem has translational invariance and the soliton/kink can therefore 
migrate. In our problem, the translational invariance is not there, due to the presence of the 
barrier and the kink solution is not free to move in space. However, the polymer on which 
the kink exists, can move, though the kink is fixed in space. Thus, the polymer goes from 
one side to the other by the motion of the kink in the reverse direction on the chain. If 
there is no free energy difference between the two sides of the barrier, then the kink moves 
by diffusion and the time of crossing 

t cross ~ A^VT^/^. If there is a free energy difference, 
then the kink moves with a non-zero velocity from the lower free energy side to the other, 
leading to t 

cross ~ N/ \/T . We also consider the translocation of hydrophilic polypeptides 
across hydrophobic pores. Biological systems accompolish this by having a hydrophobic 
signal sequence at the end that goes in ffist. Our analysis leads to the conclusion that for 
such a molecule, the configuration of the molecule in the transition state is similar to a hook, 
and this is in agreement with presently accepted view in cell biology [^]. It is also possible 
that a kink movement mechanism might operate in other biological phenomena, like protein 
folding pg 
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A. APPENDIX 
I. The eigenfunctions of the Hamiltonian 



The Hamiltonian = — ^ + cUq |l — 3sech^ (ci;o^/2)| has the following eigenfunctions 
(functions are not normalised) and eigenvalues, if n allowed to be in the range (—00, 00). 



1. Discrete States 



1) Mn) = sech(^^)3; So = -^uj^/A 

2) ^i(n) = sechl^f tanh(^^); ei = 

3) ^2{n) = {-3 + 2 cosh(a;on)} sech(^^)3; £2 = 3a;2/4 

2. Continuum States 

The continous part of the spectrum starts at uJq. The potential is refiectionless. Corre- 
sponding to an eigenvalue uJq + k'^, there are two eigenfunctions, which we write as an odd 
function and an even function. They are: 

1) ipeven{n) = 8 k + u^) cos{k u) - 3 cuq (8 P + 3uj1) sm{kn) tanh(^) 
-SOkui cos{kn) tanh(^^)2 + 15 sm{kn) tanh(^^)3 

2) ipoddin) = -8k (k^ + u^) sin(A;n) - 3a;o {8k^ + 3u^) cos{kn) tanh(5^^) 
+30 kuj^ sm{kn) tanh(^^)2 + 15 a;^ cos(A;n) tanh(^^)3 

In the limit n — > ±oo, the even function becomes like 

■061)671 (^) — ^k {Ak^ — HcUq) cos{kx) ± 6a;o {—Ak^ + <^o) sm{kx), which may be written 
as (Constant) cos{kx + S{k)), so that the phase shift 6{k) — arctan( ^ '^'^'||'^2_|_4^2| ) ■ The 
phase shift for the odd solution is just the same. Hence the total change in the density 
of states is given by An{k) = f ^ = -f (p^ + + ^^l^)- On integration, 

/q°° dkAn{k) — —3 as it should be, as there are three bound states for HK 

II. Evaluation of the Infinite Products 

1. Hairpin Crossing 



The infinite product is: 



I- hp 



k 



, n4. 



(A-1) 



where represent the eigenvalues of the continuum states of the hamiltonian H"^^ = 
(y — -§^ + a;o)and ^^.are the eigenvalues of if*, satisfying the boundary conditions at n = ±oo. 
The above product involves only the continuum eigenvalues of the two Hamiltonians. Now, 
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\nlhp = - \ J2 1^ 




where the n{k) stands for the density of states in the continuum, for the Hamiltonian if**and 
rz^p(fc) for the Hamiltonian H^. The change in the density of states is Anhp{k) = ~n{k) + 
n\p{k) and is easily evaluated from the information given in subsection [Al| . It is Anfip{k) = 



+ + sJrk^) !• Using this, we get 

15 



hp = (A-2) 



2. End Crossing 

The product that we wish to evaluate is 

1/2 

(A-3) 



I- end 



\A:>2 / 



This infinite product in the above equation is over the continuous spectra of the two Hamil- 
tonians and may be evaluated. The change in the density of states is now just half of the 
density of states density of states for the hairpin case. That is, Anf.nd{k) = ^Anhp{k). At 
the first sight, this leads to a problem, as dk Ariendik) = —3/2, instead of the expected 2 
(as has two bound states while H"^^ has none). The solution to ths is quite well known 
- if™-* has a state with eigenvalue uJq where its continuous spectrum starts, and half of this 
state is to be considered as a bound state. Then, we can write the above as: 



/ 1 f°° / 

lend = exp f -- dkAUendik) lu (uq + /c^ 



Hence we find 



lend ={f) (A-4) 
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FIG. 1. The potential energy per segment of the chain, plotted as a function of position 



FIG. 2. The potential energy along the reaction co-ordinate. The Fact is independent of the length of the chain. 
After the barrier is crossed, there is a region of width W, with W proportional to N, which is to be crossed. The 
time required to cross this regions is tcrose 



FIG. 3. The transition state 



FIG. 4. End Crossing 



FIG. 5. Hairpin Crossing 
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FIG. 6. The barrier and its inverted form. The barrier heights in the forward and backward directions are shown. 
The dotted Une represents the path that determines the activation energy 



FIG. 7. The free energy per segment of the polymer, shown as a function of the position of the segment (for the 
case where DNA is drawn through a pore). As the segment goes from the left (-ve) to right (+ve), the free energy 
changes by -AV 



FIG. 8. The full curve is the plot of the potential for the motion of the particle for < n < s, while the dotted 
curve is the potential for s < n < A'^. The particle starts at Ro at the time t = 0, moves on the full curve and reaches 
Rs at the time n = s. At this time, the potential suddenly switches to its negative. The particle then moves on this 
potential (dotted curve). The path of the particle is drawn with dashes and direction in which it moves is shown by 
the arrows. 



FIG. 9. The transition state for a hydrophilic chain with a hydrophobic signal sequence, passing through a hy- 
drophobic pore. Compare with figure 14-14 of the book by Alberts et. al. 



FIG. 10. This is not a possible transition state. 
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